function Delta = AIA_step2(I,Phi)
    F = length(I);
    Delta = zeros(1,F);
    for f = 1:F
        R = AIA_step2_R(I{f,1},Phi);
        S = AIA_step2_S(Phi);
        X = S \ R;
        Delta(f) = atan2(-X(3),X(2));
    end
    Delta = Delta - Delta(1);
end